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Abstract 

We propose an alternative evaluation of the long-distance weak annihilation (WA, also called 
one-photon exchange in this paper) contribution to the rare semileptonic 

{i = e, p) decays. This hadronic description at low energies is matched at intermediate energies 
to its short-distance counterpart in terms of quark and gluon degrees of freedom. Although the 
WA contribution does not contribute to solve the possible breaking of lepton-universality observed 
by LHCb in the —)• Ar=* * * § =(p+p“/e+e“) ratio, nor provides an important hadronic contamination 
to their decay rates, its contribution to the branching ratios (and direct CP asymmetry) of the 
B^ 7r^B£~ transitions turns out to be significant. This hadronic pollution should be taken into 
account when looking for new physics effects in decays into pions, which suggests to restrict these 
searches to squared lepton-pair invariant mass in the (1,8) GeV^ range. The interference of the 
one-photon exchange contribution with the dominant short-distance one-loop amplitude induces a 
sizable CP asymmetry in these rare decays, which calls for dedicated measurements. 
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I. INTRODUCTION 


Rare meson decays are expected to serve as harbinger for New Physics (NP) in exper¬ 
iments at the intensity frontier Better measurements or upper limits on a plethora 

of rare (semi)leptonic and radiative K and B meson decays in forthcoming experiments, 
compared to precise Standard Model (SM) predictions, eventually will provide indirect indi¬ 
cations of heavier particles with new interactions. Particularly sensitive for NP searches are 
those decays dominated by short-distance (SD) dynamics where the hadronic uncertainties 
are well under control. Conversely, precise measurements of these rare decays, combined 
with non-observation of NP effects, will furnish a better determination of flavor mixing 
parameters. 

Very recently, the LHCb collaboration has reported a dehcit in the ratio of muon to 
electron pairs produced in —)■ decays, in the (1, 6) GeV^ region for the squared 

invariant-mass of the lepton pair . This energy region is cleverly chosen as it excludes long¬ 
distance (LD) contributions associated to charmonium and light vector resonance production 
in B~ decays followed by their conversion to lepton pairs (for definiteness, hereafter we 
will focus on decays of negatively charged B mesons, unless otherwise indicated). The 
measured ratio Rk = B(B~ —)■ K~ jji~) /B{B~ —)■ K~e^e~) = (0.745to;o74 0.036) 
j^, if conhrmed in more refined measurements, would call for lepton universality violating 
interactions since in the SM = 1.0003 ± 0.0001 jsl-lsl for the energy region reported 
by LHCb j^. Possible explanations involving NP interactions have been suggested as the 
source of non-universal leptonic interactions j^. Interestingly, other anomalies have been 
reported in angular observables of related B —)■ decays lOj] which have been widely 


discussed in the literature 11 


In the SM, the semileptonic B 
SD b —)■ (s, transition jhldSl, I12M15 


—>■ P (P = K, tt) decays are dominated by the 

This elementary process is induced by the 


electromagnetic and weak penguin (Figure la), as well as W boson box (Figure lb) diagram 
contributions, which are dominated by loops involving the top quark. In these exclusive 
processes, LD chiral corrections to the hadronic matrix element have been computed for 
soft momenta of light pseudoscalars (i.e. high-momentum q of dilepton pairs), where they 


amount to rather large 0(20,30%) corrections 16|. Our point here is that, although this 


kind of corrections will be much smaller on the other energy end (large recoil region for 
the heavy-to-light B-meson form factors 0 ). the high precision of measurements as well as 
the sharp predictions of the dominating SD contributions suggests LD effects associated to 
resonance dynamics at the GeV scale be taken into account as well for these more energetic 
P mesons. Contributions of four-quark operators giving rise to tree-level WA amplitudes 
with virtual photon emission are also possible for exclusive B~ —)■ P~i~^i~ processes and 
are the main concern of this paper. 

As far as we know, the one-photon exchange Feynman diagrams (arising from WA op¬ 
erators) shown in Figures Ic-ld provide a so far neglected LD contribution that must be 
included in the calculation of these exclusive semileptonic rare decays 0. These additional 

^ We neglect a subdominant 0(Gp) neutrino-exchange amplitude. This one was checked to give a highly 
suppressed effect in —>■ 7r^£~^£ decays | 18 l| . 
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LD contributions could in principle modify the Rp value due to kinematical effects 1^. Ac 


cording to chiral effective theories of QCD 20|, for photon four-momenta below ~ 2 GeV, 
photon-hadron interactions are best formulated in terms of the relevant (pseudo-)scalar and 
(axial-)vector hadronic degrees of freedom. As it was shown long ago 2l|, the —)• 

rare semileptonic decays are dominated by the one-photon exchange contributions analogous 
to the ones of Figure Ic and 1(F the SD top-quark loop penguin and W-box contributions in 
that case are negligible small [l|. Owing to gauge-invariance, this LD contribution vanishes 


at lowest order in chiral perturbation theory but this is not the case at higher orders 2l|, |22 

According to the respective CKM mixing factors in the SD and LD amplitudes of B~ 

decays, one can envisage a correction at the percent level for P = K but a 
slightly higher effect could be expected for P = vr depending on the limits of the inte¬ 
grated observables. It is worth noticing that WA contributions were previously considered 
for B {K*/p)i+i 


[231] and B —)■ 
scheme 


25 


24j | decays in the framework of the QCD fac- 
Recently, non-local effects in P —)■ decays have been 


26| using the method developed in Ref. |27|| and previously used in Ref. 


12 


torization (QCDf 
computed in Ref. 

for B —)■ decays for low dilepton invariant masses. At large hadronic recoil (corre¬ 

sponding to low values of the squared virtual photon momentum Af^), the light-meson 
energy (Pp) is much larger than the energy scale of hadronic binding: Ep S> Pqcd- Then, 
the virtual photon exchange between the hadrons and the dilepton pair and hard gluon 
scattering can be systematically expanded in powers of 1/Pp using QCDf. However, the 
B —)■ Py amplitude suffers from nonperturbative corrections due to uu and cc intermediate 


states, which form the p, cu, J/T, ... resonances 28(| (similarly, in D —)■ Py processes dd and 


ss can be excited and the 0 meson plays a prominent role (29(|). Therefore, a clean QCDf 
prediction shall be limited to an energy range which excludes these resonant contributions, 
conservatively 2 < < 6 GeV^. 


On the contrary, the extension of these limits req uires the inclusion of near-threshold uu 
and cc {dd and ss) resonances. Let us recall Ref. [2^ where these are accounted for studying 
the b ^ s'y processes leading to a relevant effect of the charmonium resonances in the 
inclusive decay width. Also, Ref. 29(| considered resonance effects from gg excitations in P —)• 
processes, where 0 exchange basically saturates the decay width. In the analyzed 
B~ P~i~^i~ decays, cutting at a maximum g^ value well below the charmonium region 
kills the effects of the cc resonances and we have checked that the remaining contributions 
of uu resonances are negligible. There are, however, other resonance contributions that may 
be relevant in our case: the ones where the photon is attached to the initial and hnal state 
mesons in WA contributions (Fig. Ic-ld). 


In the spirit of the phenomenological Lagrangians, we would like to focus in the g^ < 2 
GeV2 region by adding the effect of the lowest-lying light-flavoured resonances and pseu¬ 
doscalar mesons as active degrees of freedom as it is done in the (Resonance) Chiral La¬ 
grangians. This is the approach that we follow to ensure the applicability of the QCDf 
analysis of the B~ P~i~^i~ decays down to threshold. Here we treat the photon-meson 
interaction by considering the exchange of resonances in the framework of two different mod¬ 
els which have shown to give a very accurate description of experimental data for g^ values 
up to 2 GeV^. This agreement shows that up to those energies, the photon does not resolve 
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the quark structure of mesons and only probes their electromagnetic structure. In order to 
avoid potential double counting of effects from WA contributions (in the QCDf approach 
and ours), we will consider our LD WA contribution only for photons with < 2 GeV^ and 
match it to the QCDf contribution for > 2 GeV^ as discussed in section HVl 

The purpose of this paper is to quantify these expectations and evaluate the contributions 
given by diagrams in Figures l(c,d) for the B~ —)■ decays for low and intermediate 

values of the invariant mass of di-lepton pair (other diagrams present in hgure 2 vanish 
due to gauge invariance, see below). We use two different approaches for the 7i~ and K~ 
electromagnetic form factors in order to estimate the systematic error of our computation. 
For these rare B decays, it is found that the one-photon exchange diagrams do not modify 
lepton universality for the energy range measured by LHCb and that they turn out to be 
very small in the P = K case, although their effects in the P = tt rates are more signihcant. 

The WA contribution that we study has a different dependence on CKM mixing elements 
than the dominant SD one and also carries a sizable strong phase, leading to possible CP 
violating effects. The dominant source of CP violation, both in our paper and in the QCDf 
approach, arises from the interference between the SD top-mediated loop and the WA di¬ 
agrams. However, the strong phase in our case is provided by the imaginary parts of the 
light resonance shapes, while in the QCDf approach it arises from the on-shell spectator u 
quark in the B meson after having emitted the photon. We hnd measurable effects both for 
F = TT, iP on the integrated CP asymmetry and encourage the LHCb and Belle-H Collabo¬ 
rations to measure such asymmetries which could be useful to validate the models employed 
in the description of exclusive b —)• (s, d)£^£~ decays at low photon virtualities. 


II. DECAY AMPLITUDE 

In the Standard Model, observables like direct CP asymmetry in B^ —)■ P^£'^£~ have 
shown to be particularly sensitive to the interference of top dominated loops and WA ampli¬ 
tudes, which carry different weak and strong phases. Both amplitudes have been evaluated 
within the framework of QCDf where the strong phase in the WA amplitude, at relatively 
low photon virtualities, arises mainly from the on-shell u spectator quark in the B meson 
after having emitted the photon. In this Section we provide an alternative evaluation of this 
WA amplitude using the Resonance Chiral Theory within a factorization approximation, 
which we expect to be valid for <2 GeV^. In a later Section we will show how to match 
this contribution to the corresponding WA amplitude calculated within QCDf, which will 
be used for larger values. Since the strong phase in our calculation stems from the pseu¬ 
doscalar electromagnetic form factors, our approach provides an independent evaluation of 
the CP asymmetry for low and intermediate values of the lepton-pair invariant masses. 

A. QCDf description 

For completeness and later use, we will consider hrst the description of these decays 
within QCDf as given in ref. j^. Let us choose the following convention for the particle’s 
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FIG. 1. Short-distance contributions to B~ —)• decays are shown in (a) (penguin) and (6) 

{W box) diagrams. The one-photon, LD, WA contributions are shown in (c) and (d). The full circle 
denotes the electromagnetic form factor of the charged pseudoscalar mesons. Other long-distance 
structure-dependent vector and axial-vector terms (shown in fig. 2) do not contribute due to gauge 


invariance 


2l| 


momenta: 

The contributions of the diagrams in Figures (la,b), including higher order QCD corrections, 
to the decay amplitude can be written, to a very good approximation as 

G pd 


M 


QCDf 


V2' 


-ytbV*pip{q^)p^P + FA{q^)h^,ld] 

TT 


( 1 ) 


where q = pp + p- is the total momentum of the lepton pair, the subindex in the CKM 
matrix element stands for D = d,s in the case of IAS"! = 0,1 transitions. The factors 
Fv ^ Ca = 4.214 and Fa = Cio = —4.312 denote the vector and axial Wilson coefficients at 
NNLO [23| corresponding to the Og = (qLlfib){i'y^i) and Ofg = (^7^75^) operators 

in the effective weak Hamiltonian for the b —)■ qi~^i~ transition. The WA contributions 
induced by four quark operators are conventionally absorbed into the effective Fv coefficient, 
by replacing its numerical value given above by ip{q^)Fv —)■ ip{q^)Fv + F'^^. For the 
F = IT case, the product of CKM matrix elements for loops involving all up-type quarks 
(VibV*^) scale as (T(A^). However, the loop function favors the top quark as the dominant 
contribution [§ owing to the hierarchy of up-type quarks. 

^ We are not neglecting the charm and up quark contributions in our numerical analysis, however. We are 
also including subleading corrections with heavy-meson form factor ratios and other Wilson coefficients 
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Althou gh w e are aware that more refined heavy-to-light meson form factors have been 
developed 12|,l3ll], we will use the classical Ball and Zwicky form factors 32(] for our estimates 
and include the corresponding uncertainty in the errors, as discussed below. The expressions 
for the g^-dependent form factors are 

. . 2^ 0.918 0.675 ^ 

1 oo r. 1 


iK{q^) = 


g7(5.32 GeV)2 
0.0541 


+ 


g2/(6.18 GeV)2 
0.2166 


+ VK{q^), 


1 - g2/(5.41 GeV)2 ' [i _ ^2/(5 41 GeV)2]2 

where the polynomials Vp{q^) as well as the needed Gegenbauer moments at the required 
energy scale can be found in Refs. 0,0. 


B. LD weak annihilation within RyT 

In this section, we would like to consider specifically light-resonance effects in the WA 
contribution to —)■ for the 7 2 GeV^ region. We assume that in this 

energy range the relevant degrees of freedom are these light mesons and resonances (and 
not the quark and gluon fields), according to the methodology of phenomenological La- 


grangians. Low-energy phenomenology supports widely this hypothesis 30|]. Applying this 


method to the process under consideration requires further hypothesis about the strong and 
electromagnetic interactions of B mesons which, we argue, play a subleading role in our 
description. 

The decay amplitude corresponding to WA is given by 


MYd = 


WA 


( 2 ) 


where denotes the effective hadronic electromagnetic current coupled to the leptonic 

current. Gonservation of the electromagnetic current demands 




(Pb + Pp)u - 


— nip 
g2 


% 


F{q^ 


( 3 ) 


where only the first term within square brackets gives a non-vanishing contribution owing 
to the condition = 0 (for the same reason, we can also replace {pp + Pp)ijl —t 2pB^). 

The factor F{q^) encodes the information about the dynamics of weak and electromagnetic 
interactions in the hadronic blob. 

Owing to gauge-invariance, most of the diagrams (shown in Figure[2]) contributing to WA, 
vanish. To illustrate this, we insert a complete set of intermediate states in the relevant 
hadron matrix element and take into account the T-ordering of the electromagnetic and 
weak interactions, which gives rise to the leading order factorizable contributions 




to eq. o, see ref. 


+(p^mfpiY^) (Y^uriB^) . 


( 4 ) 


6 




















where the sum over virtual intermediate states X and Y that are allowed by the quantum 
numbers of the weak and electromagnetic interactions, should be understood. 

Some of the Feynman diagrams that contribute to are depicted in Fig. |2l Non- 

factorizable contributions, like the ones shown in Figure [pi, are sub-leading in the 1/Nc 
expansion and will not be considered in our calculation □ . At lower photon virtualities, 
where resonance degrees of freedom may be neglected, it was shown that non-factorizable 
(local) contributions partially cancel among themselves 211] , leaving the form factors of K 


and 71 as the dominant contributions in —)■ 


( 3 ) 


D 





FIG. 2. Some of the Feynman diagrams contributing to the effective hadronic electromagnetic 
—)• vertex. Single lines stand for pseudoscalar mesons, double lines for (axial-)vector 

resonances and wavy lines for the virtual photon (due to the spin-one nature of the weak current, 
spin-zero resonance contributions are suppressed). Filled squares denote the weak/electromagnetic 
vertex, while the empty rectangles denote the WA Hamiltonian. All contributions to the —)• 

decays (including the pointlike interactions in the first line) vanish due to gauge invariance 


2ll | or are suppressed, except diagrams (i) and (j) which contribute to the electromagnetic form 


factors of pseudoscalar mesons. Odd-intrinsic parity violating vertices are also considered in the 
calculation. 


Under the above approximations, which justify the naive factorization of the hadronic 
matrix element of the weak hamiltonian, the leading order one-photon exchange (WA) am- 

^ We recall that RyT is an expansion in 1/iVc 
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plitude corresponding to Figure |2]^i,j) can be computed taking into account that 

{0\uY{^-l5)b\B-) =--f5)u\0) = i fpppf,, (5) 

and is given by (neglecting the contribution of diagram Ic -which is of order {mp/mpY 
with respect to Id- is, however, a good approximation @) 

M LD,WA — V2GpiATra)V^,V:pfBfp- ^ 


q^{m\ — rrip) 

X [Mb {Fp{q^) - l) - [Fpiq^ - l)] , (6) 

where fx denotes the decay constant of the pseudoscalar meson X according to the PDG 
33| conventions for fK,-K,B and Fx{q^) is the electromagnetic form factor of the correspond¬ 
ing meson. Note that the pure scalar QED (point meson approximation), as well as the 
structure-dependent vector and axial-vector terms, do not contribute to the amplitude as 

Thus, the one-photon exchange contribution is sensi- 


21 


required by gauge invariance 
tive to the intermediate q^ region of the pseudoscalar form factor, and eventually to their 
resonance structure. 

Due to the vector nature of the one-photon exchange contribution, its amplitude can 
be absorbed into the contribution of the Og operator in the QCDf amplitude under the 
replacement 

'Fp(g2) - r 


ip{q)Fv — ip{q)Fv + Kprup 




where 


Kp = —Stt 


2 VubV*p) 


IbIp 


VtbV*B m% - mj 


(7) 


( 8 ) 


Note that Kp ~ (9(10 x so that its influence is governed by the ratio of CKM 

factors which is ~ (9(A°) for F = n and (9(A^) for F = K. This suggests a larger effect for 
B~ —)■ transitions but a detailed analysis of the electromagnetic meson form factors 

is needed to conhrm these expectations. 


III. PSEUDOSCALAR FORM FACTORS 


For the electromagnetic form factors of the light pseudoscalar mesons (P = tt, K) we 
have considered two approaches. On the one hand we have used form factors that are ob¬ 
tained within the frame of Resonance Chiral Theory (RyT) and, on the other hand, 
the phenomenological form factors used by the BaBar Collaboration which employs the 


Counaris-Sakurai (CS) parametrization 351]. The hrst approach has the advantage of pro¬ 


viding a low-energy behaviour complying with the chiral limit of C^D 20| , which is a must 
if we want to get close to thresholds in some of our evaluations □. Alternatively, the CS 


^ The radiation emitted from lighter charged particles is expected to dominate for low values of photon 
virtualities in agreement with our result in eq. (0). 

® In particular, let us stress that this forbids contributions from Trpy or Troiy vertices for the B~ -A 'n~Ft~ 
transition (and analogous vertices for the K case). 

® On the contrary, GS parametrizations introduce spurious phases below thresholds and the decoupling of 
excited resonances at low energies is spoiled. 
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parametrizations include more excited resonances and, for this reason, are expected to give 
a closer description of data at higher energies. In the cases at hand, the bulk of the contribu¬ 
tion (even for integrated observables starting at = 1 GeV^) will be given by the 0(1020) 
(p(770)) meson exchange in the K (vr) cases. Therefore, we should expect very similar results 
for the two approaches, RxT being more reliable for observables starting near thresholds 
and the GS for the much less important higher energy range (up to 6 or 8 GeV^ depending 
on the channel). 


The minimal Lagrangian of RxT was derived in Ref. 3^ upon the requirement of chiral 
symmetry for the pseudoGoldstone helds and flavor symmetry for the light-flavored res¬ 


onances. Ref. 
counting 


230 


34l | supplements the yPT Lagrangian at lowest order (p^) in the chiral 


= V(“X + X+) , 


(9) 


with the following interaction terms between a light spin-one resonance and the pseudoGold¬ 
stone helds 34 




( 10 ) 


The real couplings Fv. Gy and Fa are unrestricted by symmetries and encode the dy¬ 
namics. One can proceed similarly for the spin-zero resonances 3^ and include terms with 
more resonance helds and excited resonances, extending also to the odd-intrinsic parity sec- 
Within the single resonance approximation, the evaluation of the P-vector form 


361,137 


tor 

factor using the PyT Lagrangian yields 


Pp(g') = 1 + 


FvG 




P2 


M2 


( 11 ) 


By demanding a Brodsky-Lepage behaviour [38| to Fp{q^) one gets the relation FyGy = 
p2^ which gives a parameter-free expression. This naive result can be improved by in¬ 
cluding the next-to-leading order ehect in the 1/Nc counting, i. e. the meson (oh-shell) 
width according to chiral constraints (the corresponding real part of the chiral loops is also 
included), and adding the ehect of higher excitations and (in the tt case) including the 
dominant isospin-breaking ehect given by the p — u mixing. 

Specihcally for F-^{q^), we have used the parametrizations in the last two Refs, in 39| 0, 
which include three isovector resonances (p(770), p(1450) and p(1700)) and the resumma¬ 
tion of hnal state interactions encoded in the chiral loop functions. These representations 
provide good quality hts to Belle data {43 for t~ decays. Additionally, we have 


^ The conventions used in the following are standard and can be checked, for instance, in 
® The numerical values that we have employed in the present analysis are those given in these references. 
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FIG. 3. RxT and GS parametrization (BaBar fit) of the electromagnetic pion form factor as a 
function of mu = compared to experimental data from BABAR 421. 


included the characteristic p(770) — a;(782) interference appearing in the neutral channel by 
multiplying the p(770) term by the factor 


9. 




( 12 ) 


with 9p^ = (—3.3 ± 0.5) ■ 10 ^ GeV^ j4l|. Other isospin breaking corrections are neglected. 


This parametrization is compared to BABAR data ^ in Figure 2, which are available for 
the energy region 2m.n- < < 3 GeV. There, the phenomenological GS parametrization 

includes an additional p-like excitation. Tiny differences between both parametrizations are 
seen in the region where the destructive interference between the p(1450) and the p(1700) 
resonances is stronger. The comparison also shows the effect of the RyT parametrization 
lacking of the p(2250) meson, which is clearly visible in the (2, 2.3) GeV region. These minor 
differences are taken into account in the hnal quoted errors. 

Analogously, we have employed two parametrizations of FV(g^) form factor. Since this 
form factor is completely dominated by the extremely narrow 0(1020) meson, we have consid¬ 
ered the RxT prediction with only one multiplet of resonances j^. On the other hand, the 
BaBar Gollaboration 4^ has reported measurements of the kaon form factor from threshold 
up to 2.5 GeV. Since |F/f(g^)p drops by 6 orders of magnitude in going from the peak to 
2.5 GeV (see Figure 3), it should be sufficient to include only the 0(1020) resonance in its 
parametrization. 
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FIG. 4. R^T and GS parametrization (BaBar fit) of the electromagnetic kaon form factor as a 
function of mu = are compared to experimental data from BABAR 441. 


As it can be observed in Figure 3, experimental data for are reasonably well 

described from threshold up to approximately 1.3 GeV. Deviations at higher energies have 
a negligible impact in the integrated observables of rare B decays. Above this energy, other 
resonance structures with very small (and alternating) peaks and dips around the single res¬ 
onance queue can be observed. These, in turn, are well described using the parameterization 
quoted in the BaBar paper ^ (which we htted to the BaBar data), which includes two 
(/), three p and three u excitations in addition to the single lightest vector meson multiplet 
included in the RxT form factor [^. 

Because of its relevance in the study of CP violating observables, we plot in Fig. [5] the 
real and imaginary parts of the light-meson electromagnetic form factors according to the 
two different descriptions that we used. While models agree nicely for the vr case, they do 
not in the K case around the p(770) peak. Since the corresponding phase shift has been 
validated by data in the case of the RxT description, we attribute this to an incorrect phase 
of the GS kaon form factor. 


Finally, the i?-meson electromagnetic form factor contribution is suppressed by a factor 
mSp/M^ in eq. ([6]) with respect to the ones of lighter mesons. In addition, a dynamical 

® Besides, in the limit of ideal a;(782) — (/)(1020) mixing the former state does not contribute to Fpiq^) and 
only p{770) and b(1020) remain, with a prominent role of the latter. 
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FIG. 5. Real and imaginary parts of the electromagnetic light-meson form factors for P = tt 
(upper half) and P = K (lower half) as a function of mu = according to the RyT and GS 
descriptions. 



form factor at infinite momentum transfer, the shape of Fsiq'^) is given just by flavour 
symmetry, yielding 


_ ^^2^ VM2-g2_,M,r,(g2) 3{M^ - 

The only way to access the electromagnetic B meson form factor at small values of the momentum transfer, 
would be by means of elastic electron-R meson scattering. Given the difficulty to produce a B meson 

target/beam, it becomes rather complicated to test this low energy behavior of the B meson form factor. 
The procedure employed is consistent with the heavy quark mass limit, as explained in the quoted refer¬ 



ence. 
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FIG. 6. Normalized lepton-pair spectrum for the LD WA contribution to B~ —>• decays 

for P = TT (upper half) and P = K (lower half). The squared lepton-pair invariant mass is taken 
from threshold up to 1.5 GeV^. 


where Tp{q‘^) can be found in Refs. 39|. Our numerical evaluations confirm that the effect 
of this form factor will be completely negligible in rates of rare B~ mesons. 

The spectra of the lepton pair, normalized to the B~ meson decay width are plotted 
in Figure [H] for low values of the lepton-pair invariant mass. We observe that above the 
muon-pair threshold both spectra are identical. Clearly, the integrated rates in the whole 
kinematical range would exhibit a trivial breaking of leptonic universality owing to the 
lower electronic threshold and its enhancement due to the 1/q^ dependence in the LD WA 
amplitude. However, this will happen both in the well-known QCDf as well as in the new 
LD WA contributions; thus only the numerical evaluation will tell if there is any additional 
measurable breaking of universality due to LD WA effects in the considered processes. 
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IV. MATCHING OF THE RxT AND QCDF DESCRIPTIONS OF THE WA CON¬ 
TRIBUTIONS 


The LD (one-photon exchange) contribution to the B~ —)■ decays that we have 

discussed in previous sections and the one evaluated in the QCDf approach originate from the 
WA four-quark operator of the effective Hamiltonian. Then it naturally arises the question 
whether there is some double-counting between both contributions or not. 

For low enough values of {q^ the QCDf expansion in 1/Ep is the best way 

to organize the perturbative series. However, when approaching the uu threshold {q^^x < 
Mj/^) light-resonance effects need to be included explicitely as active degrees of freedom in 
the action as done in the resonance chiral Lagrangians. Within the interval from threshold 
to q^T^axy expect that chiral Lagrangians give an adequate description for low values, 
while QCDf be more appropriate for the higher region. We thus assume there is an 
intermediate energy scale where both descriptions are good approximations and search for 
this matching scale where the (complex) effective Wilson coefficient Fy = smoothly 

matches both descriptions for the tt and K meson channels. We have found such matching 
scale at g^ ~ 2 GeV^, as it can be seen in Figs. [71 

Then we have consider our LD one-photon exchange diagrams as the WA contribution 
for g^ < 2 GeV^ and the SD QGDf counterpart (where the photon is radiated by the valence 
quarks in the B and P mesons) for g^ > 2 GeV^. We have used this separation scale 
between the LD and SD (QGDf) descriptions of WA contributions in the following and a 
slight variation of the matching scale has been taken into account in the error estimates. In 
the following, whenever we use QGDf it is understood that the corresponding result for WA 
is only used above 2 GeV^. 
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FIG. 7. A smooth matching between the LD and SD (QCDf) WA contributions is appreciated in 
the plot of as a function of for P = -n (upper half) and P = K (lower half). 


TABLE I. Integrated branching ratios of B~ —>• P~i~^i~ decays for P = tt (left hand side) and 
P = K (right hand side) for different ranges. We tabulate separately the QCDf , long-distance 
WA (LD) and their interference contributions for the kinematical ranges of interest. 


p- 7T-i+e- 

B- K-i+e- 


0.05 < g2 < 8 GeV^ 

1 < g2 < 8 GeV^ 

1 < g2 < 6 GeV^ 

LD 

(9.06 ±0.15) • 10“^ 

(4.74 ±0.05) • 10-^9 

(1.70 ±0.21) • 10-9 

interf. 

(-2.57 ±0.13) • 10-9 

■ 10-^9 

(-6 ±2) • 10-1^ 

QCDf 

(9.57il;^f) • 10-9 

(8.43l(,;9i) . 10-9 

(1.90l[j;^?) X 10-^ 

Total 

(1.6ll9;if) • 10-8 

(8.69lJ;9i) . 10-9 

(1.92l9;9?) X 10-^ 


V. BRANCHING FRACTIONS FOR R 


p-l+t 


Searches for New Physics at large hadronic recoil in the B~ decays are 

restricted to the (1, 6) GeV^ range of j^. Thus we will stick to this region for this 
channel. The B~ —)■ decays have just been observed [46| and such studies have not 

taken place yet. We therefore include our results both for a range starting at = 0.05 
GeV^ (basically the muon threshold) and at 1 GeV^. In either case we cut the phase-space 
integration at 8 GeV^ to avoid the charmonium region. In Table [I] we show the corresponding 
branching ratios of LD WA, QGDf and their interference contributions by considering these 
different kinematical integration domains. Values for required input parameters are taken 
from the PDG 3^ except for the GKM matrix elements which come from Ref. 

Several comments concerning these results are in order: 
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Our results for the QGDf contribution to the B tt branching ratio are higher 


than those in ref. 15| because of the different heavy-meson form factors employed. 


Specihcally, in this analysis form factors parameters were htted to reproduce B~^ 


15 
























data, resulting in smaller QCDf contributions than in other analyses 13|, [14 


or ours. 


Another source of difference in the QCDf contributions comes from the updated inputs 
while older PDG values were employed in earlier analyses. As a 
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we are using 

result, our numbers for the vr case are larger by 

We have not performed a dedicated study of the errors of the QCDf contributions. 
Errors quoted in table |T] are obtained rescaling the errors in Refs. 151 and according 
to the different central values obtained by them and us. As discussed extensively in 


these references (see also 12|), the dominating error for the K case comes from the 


heavy-meson vector form factor, while in the tt case the choice of the renormalization 


scale fib basically saturates the overall uncertainty (see, however, jll| ) 


Our study of the LD WA contributions to B~ —)■ P~£^i~ has been performed with 
two different sets of form factors in the (1, 6) {K) and (1, 8) (vr) GeV^ ranges. In the 
above-GeV intervals, the error has been estimated from the difference between these 
predictions. When including the region immediately above threshold we have only 
employed the set of chiral-based form factors estimating the error as the difference 
between the results obtained using dispersive form factors and a Guerrero-Pich-like 


39|. Analogous procedure has been employed in order to obtain our 


resummation 
results for GP violation in the next section. 


• The violations of lepton universality induced by kinematical effects on Rk and Rtt 
are always given by the QGDf contribution. The LD WA modihcation is -in all the 
considered energy ranges- smaller than the error of the QGDf contribution. Therefore, 
Rk = 1.0003(1) in the (1,6) GeV^ range and R^^ = 1.0006(1) in the (1,8) GeV^ 
range. 


When we add the contributions of values above 8 GeV^, ref. 15| , our branching fraction 
corresponding to the full kinematical domain becomes B{B~ tt~£^£~) = (2.61q 3) ^ 10“^ 
These results can be compared to available data from Refs. 
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B{B- 7r-/i+/i") = (2.3 ± 0.6 ± 0.1) x 10"®, (14) 

B{B- R-e+e-) = (1.56l|];?|?) x 10"^ for 1 < < 6 GeV^. (15) 

Our results and experimental data agree within error bars0. In the B~ —)■ decays, 

current errors on the (completely dominating) SD contribution do not allow to tell whether 
there is a tension between SM prediction and the LHGb measurement or not. With the 
smaller error expected on the branching fraction of the B~ —)■ 7i~Ri~ decays from the next 
run of LHG measurements, one might be able to notice a tension between SM predictions 
and data. It must be noted that a reduction of the current error to less than a half will be 
able to pinpoint the LD WA contribution to these decays that we have been discussing. 


Recent results reported in Ref. 


47j| give a smaller BR for the muon case, but still consistent with the 


previous measurement. 
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VI. CP VIOLATION 


A direct CP asymmetry can be generated because the LD one-photon exchange and the 
QCDf contributions to the amplitude of the —)■ P^£^i~ decays have different weak and 

strong phases, as it can be easily checked from eqs. ([T]) and (jb]) 0 More precisely, the strong 
phase required to generate CP asymmetry in our case, arises from the LD WA contribution. 
It will dominate at low photon virtualities, since it stems from the hnite decay widths 
of vector mesons that describe the electromagnetic structure of charged mesons. Thus, a 
(partially integrated over a hnite range) CP asymmetry 0 


^cp{P) — 


T{B+ P+i+i-) - T{B- p-i+i-) 


(16) 


T{B+ P+i+i-) + T{B- P-i+i-) ’ 

can be generated from the interference of diagrams shown in Figure 1. By inserting the 
amplitudes of LD WA and QCDf contributions in the previous expressions, it can be shown 
that the width difference has the form (as seen in table [H interferences are much smaller 
than the QCDf and LD WA contributions): 

Acp = r{B+ p+i+r) - r{B- p-i+r) 

= -32a^GlfpfBlm{VtbV;j,V:,VuD} 

1 


X / dq / ds 


12 ‘ 


2{Pb ■ P+){Pb ■ P-)- 


m. 




2 

1 )]} 


(17) 


‘ — rrip) 

xlm {U(l")Fv{q^) [M| {Fp{q^) - l) 
where Si 2 = {pk +P+)^- 

According to the matching of the LD and QCDf descriptions of the WA contributions, 
the latter will also violate CP for q^ > 2 GeV^, arising from on-shell radiating light-quarks. 
Finally, there will also be another source of QCDf CP violation from light qq in loops 24 . 
Taking all of them into account it leads to: 

r (16.1 ± 1.9)%, for P = TT, 0.05 < q^ < 8 GeV^ 

Acp{P) = < (7.8 ± 2.9)%, for P = TT, 1 < < 8 GeV^ (18) 

[ (-1.0 ± 0.3)%, for P = a:, 1 < g2 < 6 GeV^. 

In the hrst energy range for the vr meson case, 83% of Acp{ti) has a LD WA origin, while 
this reduces to 31% in the second. In the case of kaons, for the 0.05 < q^ < Q GeV^ 
range, the amount of GP violation basically doubles but is completely dominated by LD 
WA contributions (which is already 70% in the 1 < < 6 GeV^ interval). The quoted error 

in our results stems from the systematic error attributed to parametrizations of the light 
meson electromagnetic form factors. 

In ref. 2^, using QGDf, the Acpiir) asymmetry was predicted to be larger than ours (for 
comparison, see Tab le ITT] for three ranges of the lepton-pair invariant mass), while the results 
in the recent paper [26| lie somehow in between of both predictions. Note that in the QGDf 


A possible large CP violation in the —>■ decays induced by WA one-photon exchange contri¬ 

bution was proposed for the first time in ref. 2^ within the framework of QCDf. 

Other CP violating observables can be analyzed analogously. 
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TABLE II. Our results for Acp{t^) (in %) are compared to those in ref. 24 f or different energy 
ranges. We note that in the (1,6) GeV^ range the recent result in Ref. 26|, is in 

agreement with both determinations, although closer to ref. 24 ] . 


i^mini Qmax) 

Ref. [24] 

Our results 

(1,8) GeV^ 

(1.6) GeV^ 

(2.6) GeV^ 

13 ±2 

16 ±2 
1311 

7.8 ±2.9 

9.2 ± 1.7 

7.7 ±0.5 


approach the strong phase for low values of is dominated by the photon emission off the 
spectator u quark (in the B meson) in the WA diagram, which may be on-shell after having 
emitted the photon (for larger values of a smaller strong CP phase arises from light qq 
in loops from form factor B ^ tt contributions). We note that this CP violating effect is 
maximal towards the q^ threshold, where the applicability of QCDf is more questionable. 
Measurements of binned CP asymmetries in this and the —?• {V = p, K*) decays 

may be sensitive to the modelization of the WA contributions at long distances. 


We note that the signs of the results in refs. 2J, |2^ have been switched in Table ITT] 


because our convention for dehning Acp{t^) is opposite. Within our approach, the solely 
contribution of the B-meson electromagnetic form factor to the CP asymmetry is completely 
negligible, at the 0(10“^) level. 

This (P(few %) CP violating hgures shall enhance the case for their measurements. Cur¬ 


rent values at the PDG are well compatible with zero 33 


in B^ —)■ (and also in 

—)■ ) decays, while this observable is not yet reported in the B^ —)■ 

case. The most recent measurement reported by the LHCb collaboration for the integrated 
CP asymmetry, namely Acp = 0.11±0.12±0.01 47|, is still consistent with the different 
theoretical predictions. 

Therefore, despite the fact that the tree-level one-photon exchange diagrams give a small 
contribution to the decay rates (specially for the K case), they can generate a non-negligible 
CP asymmetry within the Standard Model. This CP asymmetry altogether with measure¬ 
ments of the decay rates can be used as a test of New Physics in the rare —)■ 

decays. This makes us emphasize the need of a dedicated measurement of these observables 
in the next LHC run at LHCb and in the forthcoming Belle-H experiment. 


VII. CONCLUSIONS 

In this paper we have considered the one-photon exchange contribution to the rare 

decays, with P = n or K. Its effects in the decay rates of the P = K case turn 
out to be of order 1% with respect to the (top quark loop dominated) SD contribution for 
the range 1 < < 6 GeV^ of the squared lepton pair invariant mass. We do not foresee 

forthcoming measurements being sensitive to this contribution in the near future. On the 
contrary, this fact conhrms the suitability of this range for new physics searches. 

In the case of a vr^ meson in the hnal state, the corresponding effect turns out to be 
signihcant in integrated observables starting close to threshold. This suggests to take -in 
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analogy to the case of a final state with K - the range 1 < < 8 GeV^ for precision 

measurements, since the LD WA contribution is reduced to less than 10% with a negligible 
uncertainty in that interval. On the other hand, more rehned measurements of the fully 
integrated branching fraction for this decay could be sensitive to our contribution once the 
error is reduced below a half of the current uncertainty. 

Interestingly, the different weak and strong phases of the QCDf and LD WA (one-photon 
exchange) contributions are capable to generate a CP asymmetry. Again, this CP asymmetry 
is large in the case of a pion in the hnal state for 0.05 < q‘^ < 8 GeV^ values of the lepton-pair 
invariant mass, but also sizable and worth to measure in the 1 < < 8 GeV^ interval. For 

the kaon case, the range 1 < 5 '^ < 6 GeV^ is optimal for such a search. Our CP violation 
results are smaller than those obtained within QCDf because of the different description of 
the WA amplitudes at low energies. Future measurements shall be sensitive to this kind 
of contribution and shed light on its appropriate description. More rehned measurements 
of this CP asymmetry and of the magnitudes of the decay rates at LHCb and future B- 
superfactories can provide another non-trivial test of the Standard Model or may furnish 
indications of New Physics. 
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